Composite domain walls in flat nanomagnets: the magnetostatic limit 
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We discuss the structure of the so-called "vortex" domain walls in soft magnetic nanoparticles. A 
wall of this kind is a composite object consisting of three elementary topological defects: two edge 
defects with winding numbers —1/2 and a vortex with a winding number +1 between them. We 
provide a qualitative model accounting for the energetics of such a domain wall. 



Magnetic nanoparticles generate considerable interest 
as prospective building blocks for nonvolatile random- 
access memory ^1.]. Storing of information bits is made 
possible by the existence of two (or more) stable magnetic 
configurations. Switching between the two stable states 
can be achieved by applying an external magnetic field 
or by injecting current. In either case the switching pro- 
cess proceeds through the formation of complex transient 
patterns of magnetization 2] . Building fast and reliable 
magnetic memory thus requires a thorough understand- 
ing of magnetization dynamics in these nanomagnets. 

In magnetic nanoparticles with the geometry of strips 
and rings the switching creates domains in which the 
magnetization is forced by the magnetostatic forces to 
be parallel to the edge. These domains grow and shrink 
at the expense of one another until one of them occupies 
the entire sample This process can also be viewed 

as the creation, propagation, and annihilation of domain 
walls. Thus the question of local stability and dynamics 
of the magnetic configurations can be answered by study- 
ing the static and dynamic properties of domain walls. 

Domain walls in submicron rings and strips have a con- 
siderably complex structure. For instance, McMichael 
and Donahue |j| have observed, among others, configura- 
tions termed "transverse" and "vortex" domain walls. In 
a previous paper Q we pointed out that the transverse 
walls are composite objects built from two elementary 
topological defects. In the limit where exchange inter- 
action is the dominant force (thin and narrow strips or 
rings) , the two elementary defects are vortices with frac- 
tional winding numbers -1-1/2 and —1/2; these defects are 
confined to the edges because of their fractional topolog- 
ical charges. To our knowledge, these edge defects were 
first discussed by Moser |^ and Kurzke 0. 

In this paper we analyze the structure of domain walls 
in a strip in a different limit where the dominant forces 
are magnetostatic. This limit, in which both the thick- 
ness and width of a strip exceed the exchange length 
A = Aj i-ioMn , is relevant to the ongoing experimen- 
tal studes 2j j| . The nonlocal nature of the dipolar inter- 
actions makes the analysis considerably more difficult. 
Valuable information concerning the global structure of 
a domain wall is provided by topological considerations 
Under very general circumstances, a domain wall 
in a nanostrip is a composite object containing several 
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FIG. 1: Top: a magnetization configuration free of bulk mag- 
netic charges, —V • M = 0, and containing two —1/2 edge 
defects and a -f-l vortex in the middle. Parabolic segments 
of Neel walls are shown by dashed lines. Bottom: a head- 
to-head vortex wall obtained in a micromagnetic simulation 
using OOMMF [TH in a permalloy strip of width w — 500 rim 
and thickness t — 20 ran. 



elementary topological defects, some of which reside in 
the bulk and others at the edge. The topology restricts 
possible compositions of a domain wall thus providing a 
basis for selecting appropriate trial states. 

In a companion paper 9] we have identified three ele- 
mentary topological defects that survive the transition 
from the exchange limit to the magnetostatic regime: 
the vortex (winding number n — -1-1), the antivortex 
{n = —1) and one of the edge defects {n = —1/2). Bare- 
bones versions of these defects can be constructed using 
van den Berg's method wherein the exchange energy 
is initially neglected and the magnetostatic energy is min- 
imized absolutely by preventing the appearance of mag- 
netic charge — V ■ M. (The meat can be grown by includ- 
ing the exchange interaction perturbatively.) Only the 
vortex retains its original shape; the antivortex morphs 
into a cross tie (two intersecting 90-degree Neel walls); 
the —1/2 defect looks like a cross tie pinned at the edge. 
The -1-1/2 edge defect likely has a high magnetostatic en- 
ergy. We have also shown that the simplest domain 
wall in this limit is expected to contain two —1/2 edge 
defects and a vortex. We next discuss the structure of a 
domain wall in this limit. 
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FIG. 2: Top: a model vortex wall with the vortex absorbed by 
the edge. Bottom: a transverse wall observed in a numerical 
simulation (permalloy, w — 80 nm, t — 20 ran.) 



FIG. 3: Top: a model vortex wall with the vortex near the 
edge. Bottom: a similar configuration observed in a numerical 
simulation (permalloy, w — 200 nm, t = 20 nm.) 



A head-to-head domain wall carries a fixed nonzero 
amount of magnetic charge {2MQtw in a strip of width 
w), with a finite density in the bulk, — V-M 0, or at the 
film edge, n-M (or, most likely, both). Therefore van 
den Berg's method is not, strictly speaking, appficable. 
Nonetheless, an examination of the detailed structure of 
a vortex wall (bottom panel of Fig. ^ p| shows that it 
indeed contains two —1/2 edge defects and a -1-1 vortex 
in the middle. The edge defects share one of their Neel 
walls; the vortex resides at its middle point. 

In what follows we consider a model of the vortex do- 
main wall that is free of bulk magnetic charge. Thus all 
of the charge 2MQtw is expelled to the edges. Under this 
restriction it is possible to construct a vortex domain wall 
by piecing together the vortex and two —1/2 edge defects 
as shown in Fig. Q The resulting structure contains do- 
mains with uniform and curling magnetization. In a strip 
\y\ < w/2 with the shared Neel wall x = y and the vor- 
tex core at (u, v), the two curling domains in the regions 
±v < ±y < w/2 are separated by parabolic Neel walls 
(x — w)^ = (2y ± w){2v ± w) from domains with horizon- 
tal magnetization; they also merge seamlessly with other 
uniform domains along the lines x = v and y — v. 

This trial state exaggerates the accumulation of sur- 
face magnetic charge thereby overestimating the normal 
component of magnetization n • M at the edge. Nonethe- 
less it captures the major features of a vortex domain 
wall. 

The construction of a domain wall out of the three 
defects is not unique and has (at least) one degree of 
freedom: the vortex can be placed anywhere along the 
shared Neel wall x — y. When the vortex core reaches 
the edge, it is absorbed by the —1/2 edge defect. Their 
fusion creates an edge defect with the winding number 
+1 — 1/2 = +1/2 (Fig.EJ). As can be seen from the figure, 
the +1/2 defect is rather extended (length 2w) and, in 
accordance with our previous remarks, contains magnetic 



charge (all of it at the edge in this model). This structure 
is topologically equivalent to the transverse domain wall 
in the exchange limit 5], where both edge defects are 
pointlike. The transverse walls observed in experiments 
3] and simulations Q are midway between these two 
extremes: the +1/2 defect definitely has a wider core, 
although its extent is less than 2w. 

To determine the equilibrium configuration of the com- 
posite wall we computed the total energy of the compos- 
ite wall and minimized it with respect to the vertical 
coordinate v of the vortex. The energy is the sum of the 
following terms. 

The magnetostatic energy coming from the Coulomb- 
like interaction of the magnetic charges spread along 
the edges with the line densities Ai^2 = tM-ni^2 = 
±tMQsm9 for the upper and lower edges, respectively. 
It includes the interaction of magnetic charges on the 
same edge and on different edges: 



^ r At(x) Xjjx') 
8tt J \x~x'\ 



dx dx' , 



= f ^M^M^dxdx', (la) 

The total magnetostatic energy -^y of the 

order Aw{t^ /X^)\og{w/t). It represents the dominant 
contribution in sufficiently wide and thick strips. 

The energy of the Neel walls can be computed as a line 
integral 



£^wans = t ai£) di, 



(lb) 



where d£ is a line element of the wall. The wall surface 
tension a depends on the angle of rotation across the 
wall, which stays at 90 degrees along strmght segments 
and varies along parabolic ones. See Ref. 9 for details of 
the calculation. This term is of the order Atw/X. 
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In addition to these two wall configurations, which 
have been previously discussed in the literature, we have 
found two metastable states that corrsepond to local min- 
ima of energy E(v). One of them has the vortex core 
rather close to but not exactly at the edge (top panel 
of Fig. 13) . We have observed domain walls of this kind 
in numerical simulations (bottom panel of Fig. The 
other metastable state occurs when the energy curve 
E{v) has two symmetric minima around v — 0. This 
is a vortex wall with the vortex core slightly off center. 
Because the off-center minima are rather shallow it may 
be difhcult to observe such states in practice: even slight 
imperfections of the nanoparticle can change the poten- 
tial landscape E{v). 



FIG. 4: Energy of the vortex domain wall as a function of the 
vortex position u at a fixed strip width u) = 50 nm for several 
thicknesses t. 

The exchange energy 

Ee^che^nsc = At / {VOfdxdy, (Ic) 

Jn 

where fl is the area around the vortex where magnetiza- 
tion curies. This term is of the order Atlog (w/X). 

An investigation of the phase diagram and a quanti- 
tative comparison of this model with the numerical and 
experimental results is currently under way. In what fol- 
lows we report some preliminary findings. 

The evolution of the energy curve E{v) at a fixed width 
w and varying thickness t is shown in Fig.^ For substan- 
tially wide and thick strips, the one and only minimum 
of energy is achieved with the vortex in the middle of the 
strip, in agreement with numerical simulations 0. 

As the cross section decreases, a local minimum devel- 
ops with the vortex core at the edge of the strip. In this 
configuration (Figl^J the vortex and the —1/2 edge de- 
fect have merged to form an extended +1/2 edge defect. 
The configuration is highly reminiscent of the transverse 
wall Q that is known to coexist with the vortex wall 
over a range of cross sections • The transverse wall be- 
comes a global minimum of energy when the cross section 
becomes small enough. The vortex wall remains locally 
stable until it becomes a local maximum of energy. 

(The reader should note that the curve for thickness 
t = 1 nm shown in Fig. ^ is only an extrapolation: the 
energetics of the Neel wall in very thin films is a nonlocal 
problem Q to which our estimate based on Eq. (4) in 
the companion paper does not apply. Nonetheless, 
the overall trend refiected by the shape of the curve is 
correct: in the exchange limit domain walls containing 3 
defects are locally unstable.) 



The simple model of a domain wall in the magneto- 
static limit presented in this paper shows a qualitative 
agreement with observations. Its quantitative compar- 
ison with available experimental and numerical data is 
in progress. At a minimum, the model provides an in- 
sight into the nature of the vortex domain walls in a 
regime relevant to experiments. It corroborates our ear- 
lier suggestion 5;] that domain walls in nanostrips can 
be viewed as composite objects, which may be helpful in 
understanding their dynamical properties. 
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